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Abstract 

It  has  been  shown  possible  experimentally  to  produce  nulls  in  the 
pattern  of  prime  focus  reflector  antenna  using  discs  mounted  on  the  dish. 

A  model  of  this  antenna  system  is  developed  to  evaluate  optimal  configura¬ 
tions  and  ideal  performance. 

Aperture  integration  is  the  method  of  analysis  used.  Discs'  effects 
are  modeled  as  a  phase  shift  on  the  aperture.  No  secondary  effects,  such 
as  diffraction,  are  considered. 

Based  on  the  model  developed,  guidelines  are  presented  for  antenna 
design.  Computer  code  was  written  to  implement  the  model  and  a  prediction 
of  an  antenna  system's  performance  is  presented. 


NULL  GENERATION  USING  DISCS  ON  A  REFLECTOR 


Introduction 


Historical  Background 

An  adaptive  antenna  system  may  be  viewed  as  a  spatial  filter  that 
provides  an  improved  signal  to  noise  ratio  (the  noise  may  be  any  unwanted 
signal)  compared  to  a  conventional  antenna.  The  signal  to  noise  ratio 
may  be  improved  by  either  increasing  gain  in  the  direction  of  the  desired 
signal,  or  decreasing  gain  (cancellation  or  nulling)  in  the  direction 
of  the  noise. 

Adaptive  antenna  systems  have  historically  been  arrays  whose  outputs 
are  processed  by  algorithms  capable  of  steering  "main  beams,"  nulls, 
or  both.  Modern  arrays  provide  excellent  performance  in  a  multitude 
of  applications  due  to  the  inherent  flexibility  of  a  software-driven 
system. 

Arrays,  however,  have  characteristics  such  as  high  complexity, 
weight,  cost,  low  maximum  gain,  etc.,  which  may  render  them  unsuitable 
in  some  situations.  Reflector  antennas  alternately  are  of  simpler 
construction  and  offer  high  maximum  gain  per  unit  weight.  Therefore, 
a  method  of  controlling  the  far  field  pattern  of  a  large  reflector  to 
produce  improved  signal  to  noise  ratio  (by  effectively  producing  nulls 
in  the  direction  of  interference  sources)  appears  to  offer  a  new  flexibility 
to  a  highly  efficient,  and  widely  used,  antenna  system. 

This  work  was  originated  by  Dr.  J.  Leon  Poirier  and  Daniel  Jacavanco 
at  the  Electromagnetic  Sciences  Division,  Rome  Air  Development  Center, 


Hanscom  AFB,  Massachusetts.  Extrapolating  from  work  on  large,  flexible 
reflectors  controlled  electrostatically  to  achieve  a  perfect  shape 
(i.e.,  parab^’ic,  etc.);  they  theorized  that  by  distortion  of  the 
reflector  in  the  same  manner,  pattern  nulling  may  be  possible.  Work  by 
Havens  (Havens,  1983)  at  the  Air  Force  Institute  of  Technology  showed 
this  to  be  theoretically  possible.  However,  the  control  problem  of  a 
perfect  shape,  let  alone  selective  distortion,  appears  to  be  too  burdensome 
for  practical  use  (Lang,  1982:666). 

In  the  summer  of  1983  Jacavanco  demonstrated  it  was  possible  to 
produce  nulls  in  the  pattern  of  a  rigid  reflector  by  mounting  discs  on 
the  reflector  dish  and  adjusting  their  distance  from  the  dish  (Jacavanco, 
1984c).  This  method  of  reflector  distortion  appears  to  offer  an 
inexpensive  and  simply  constructed  antenna  system  capable  of  null  generation 

While  the  technique  had  been  experimentally  demo  .  trated,  there  was 
little  theoretical  basis  for  the  mechanisms  producing  the  nulls.  There¬ 
fore,  this  study  was  undertaken  to  attempt  to  develop  theoretical  back¬ 
ground.  Ideally  this  work  will  produce  a  means  of  calculating  optimal 
configurations,  bandwidths,  etc. 

Approach 

Aperture  integration  was  selected  as  the  most  straightforward  tech¬ 
nique  of  analyzing  the  antenna  system  with  discs  mounted  on  the  dish. 

While  it  is  well-known  that  this  technique  is  limited  in  effectiveness 
to  the  "vicinity"  of  the  main  beam,  it  is  in  this  area  that  null 


synthesis  has  been  shown  effective  and  is  of  interest. 


Coordinate  System 


In  order  to  avoid  confusion  a  standard  coordinate  system  is  used 
throughout  this  paper.  The  focal  point  of  the  parabolic  reflector  is 
defined  as  the  origin  with  z=0  the  focal  plane  and  the  antenna  located 
in  the  z<0  half-space.  Spherical  coordinates  are  used  interchangeably 
with  rectangular  where  they  are  more  convenient.  Coordinates  that  refer 
to  the  aperture  plane  (source)  will  be  primed  any  time  there  is  any 
ambiguity. 

Electric  Fields 

Equations  for  the  far  field  electric  field  generated  by  an  aperture 
as  given  by  Stutzman  and  Thiele  (Stutzman  and  Thiele,  1983:383)  are 

jBr 

Ee  "  je  fer  (Px  C0S  ♦  +  Py  sin 
ejer 

E,  =  jB  ~ -  cos  9  (P  cos  <J>  -  P  sin  (J>) 

<t>  2-rrr  y  x 

Since  relative,  not  absolute,  field  strengths  are  of  interest  (relative 
to  main  beam  gain),  Eqs  (1)  and  (2)  are  simplified: 

En  «  P  cos  <j>  +  P  sin  $  (3) 

9  x  y  v  ' 

E,  «  cos  9  (P  cos  <J>  -  P  sin  <J>)  (4) 


Polarization 


The  choice  of  polarization  convention  is  not  straightforward.  As 
discussed  by  Ludwig  (Ludwig,  1973),  no  one  convention  is  universally 
accepted  or  optimum  for  all  situations.  The  spherical  coordinates  for 
the  far  field  (as  in  Eqs  (3)  and  (4)),  however,  cause  serious 
difficulties  when  attempting  to  speak  practically  of  co-  and  cross¬ 
polarization  in  the  far  field  (linear  polarization  is  assumed  throughout 
this  paper).  Therefore,  a  simple  approach  will  be  used  in  order  to  avoid 
confusion. 

Both  source  and  far  field  patterns  will  be  expressed  in  rectangular 
coordinates.  The  aperture  will  be  defined  in  the  z=0  plane  and  the 
far  field  pattern  broken  down  into  x  and  y  components.  Admittedly, 
in  the  far  field  the  fields  will  have  z  components,  however  these  will 
be  small  in  the  area  of  interest  and  will  not  be  computed.  This 
component  is  readily  available  from  the  requirement  that  the  total  field 
be  orthogonal  to  the  spherical  radial  unit  vector. 

This  choice  of  reference  polarization  has  an  additional  advantage 
in  that  it  simplifies  Eqs  (3)  and  (4).  The  far  field  electric  field 
components  are  now  functions  of  only  the  like  component  of  the 
P  vector. 


I  Ei  |  2  *  cos2  0  {(RJPjD2  +  (W)2)  i  -  x,y  ( 

(Note  that  the  electric  field  is  squared  to  provide  a  power  measure.) 
P-Vector  Evaluation 

The  P-vector  used  in  Eqs  (1)  and  (2)  is  generated  by  integration 
of  the  electric  field  over  the  aperture.  The  formulation  as  given  by 


Stutzman  and  Thiele  (Stutzman  and  Thiele,  1981:426)  for  a  circular 

aperture  is 

•2n  ra  _  _ 

E  (r')  exp  {j8r^  sin  6  cos  ($  -  40)  r'dr'dip'  (6 

Jo  J0  a 

where  Ea(r^)  is  the  electric  field  on  the  aperture  of  radius  a,  (9,<f>)  are 
far  field  coordinates,  and  are  source  (aperture)  coordinates. 

The  double  integral  in  Eq  (6)  presents  a  problem  since  the  discs 
produce  discontinuities  in  the  phase  of  the  aperture  electric  field. 

The  discontinuities  prompt  the  use  of  numerical  techniques  to  evaluate 
the  integral.  The  method  is  taken  from  Abramowitz  (Abramowitz  and  Stegun, 
1972:894) 

n 

//f(x,y)  dxdy  =  v  Z  w  f(x. ,y.)  +  R  (7 

i=l  1  11 

where  v  is  the  area  of  the  sections  of  the  partitioned  surface,  the  w^'s 
are  weighting  functions  given  for  a  particular  integration~scheme,  and 
R  is  the  remainder  term  which  will  be  neglected.  Again,  since  the  v 
is  a  constant  multiplier  on  all  terms  and  only  relative  numbers  are  of 
interest,  it  is  dropped,  yielding 

n 

P  «  Z  w.  E  (r.)  exp  (jgr  sin  9  cos  (♦-$.)}  (8 

i-1  1  3  1  1 


Letting 

Vri}  =  ejYi 


(9 


Euler's  formula  is  used  to  break  P  into  real  and  imaginary  parts  yielding 


(10) 


Re[P]  «  Z  w.  |E  (r  )  ]  cos  {jf3r.  sin  0  cos  (<f>  -  $  )  +  y , } 
.  .  1  3  1  1  JL  3. 

1*1 


I 


n  __ 

Im[P]  «  E  wi  |Ea(r1)|  sin  {jSri  sin  9  cos  ($  -  4^)  +  (11) 


Eqs  (10)  and  (11)  are  the  basis  for  the  integration  scheme  as  implemented. 
Four  variables  of  summation  were  used — the  real  and  imaginary  parts  of 
P^  and  the  real  and  imaginary  parts  of  P^ — to  sum  over  the  aperture  plane 
for  a  given  9  and  $  (far  field).  These  quantities,  when  used  in  Eq  (5), 
yield  a  quantity  proportional  to  the  magnitude  of  the  electric  field 
squared,  which  is  normalized  to  the  value  obtained  for  9=<p=0,  yielding 
a  far  field  power  pattern  point. 


Antenna  Model 


The  antenna  model  is  a  prime-focus  parabolic  reflector  fed  by  a 
point  source  producing  a  parabolically  tapered  aperture  field.  The 
polarization  of  the  feed  is  assumed  to  be  linear  with  the  only  cross- 
polar  components  of  the  aperture  electric  field  those  induced  by  the 
reflector  dish.  The  feed  is  assumed  to  block  a  circular  area  in  the 
middle  of  the  aperture — essentially  setting  the  electric  field  to  zero. 

No  direct  radiation  from  the  feed  is  considered. 

The  disc's  effect  on  the  aperture  field  is  modeled  as  a  phase  shift 
across  a  circular  projection  of  the  disc  on  the  aperture  plane  (see  Figure  1) 
Amplitude  variation  in  the  field  due  to  the  different  path  lengths  is 
not  taken  into  account.  The  phase  shift  of  the  aperture  field  is  estimated 
two-dimens ionally,  as  shown  in  Figure  2,  by  computing  the  distance  from 
F-D^-A^,  and  subtracting  it  from  twice  the  focal  distance  of  the  dish 
(the  normal  path  length  from  the  feed  to  the  aperture  plane) . 

Two  approximations  that  are  made  here  are  worth  mentioning.  First, 
the  ray  tracing  approach  used  does  not  take  into  account  the  actual  geo¬ 
metrical  optics  path  from  the  disc  to  the  aperture  plane.  All  rays  reflected 
off  the  disc  are  assumed  to  travel  on  a  path  normal  to  the  aperture  plane — 
thus,  the  spreading  from  the  disc  is  not  accounted  for. 

Second,  the  disc  model  is  two-dimensional  and  does  not  take  into 
account  the  curvature  of  the  disc  in  the  third  dimension.  Thus,  the 
phase  shift  has  minor  errors  except  on  the  radial  line  running  through 
the  center  of  the  disc.  In  both  cases  it  is  assumed  the  errors  introduced 
are  minor  for  discs  which  are  relatively  small  in  relation  to  the  disc. 
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Figure  2.  Disc  Model 


Results 


This  method  of  breaking  down  the  aperture  integral  produces  simplified 
equations  that  allow  insight  into  the  effects  produced  by  the  discs  and 
methods  to  control  these  effects.  Eq  (5)  shows  that  the  co-polar  pattern 
depends  only  on  the  P  vector  component  of  the  same  polarization.  Thus, 
only  two  of  the  four  variables  summed  over  the  aperture  plane  control 
the  co-polarized  pattern  strength.  The  two  variables  of  interest  are 

Re[Pj]  =  Z  wt  |  Ej  (r^)  |  cos  [  Sr ±  sin  0  cos  (<t>  -  4>±)  +  Y±]  (12) 

Im[P  ]  =  E  wt  |E  (r^l  sin  [6^  sin  6  cos  (<f>  -  4^)  +  y^  (13) 


where  j  is  the  co-polarized  unit  vector  and  i  references  the  points  being 
summed.  With  no  discs  present  let  y^=0  for  all  i  and  examine  the  magnitude 
of  the  P^'s.  The  magnitude  of  the  aperture  electric  field  is  assumed 
to  be  a  function  only  of  the  radial  position  on  the  aperture  plane.  Letting 
radial  position  on  the  aperture  plane  (r^)  equal  a  variable  number  of 
wavelengths,  k,  and  solving  for  the  maxima,  minima,  and  zeroes  of  the 
modulation  produces 


k 


_ n _ 

4  sin  9  cos  (4>  -  4>^) 


n  =  0,1,2,... 


(14) 


For  a  four-foot  aperture  at  a  frequency  of  2.8  GHz,  a  plot  of  the 
lines  of  the  maxima,  minima  and  zeroes  for  the  (a)  cosine  and  (b)  sine 
terms  is  shown  for  <j>=0  and  9=10°  in  Figure  3.  Figures  4  and  5  show  a  plot 
of  the  modulation  (an  electric  field  of  unit  amplitude  uniformly  over 


the  aperture  is  assumed  with  a  six-inch  feed  blockage)  taken  along  the 
x-axis  of  the  aperture.  This  shows  that  the  imaginary  part  of  the  P^'s 
components  are  odd-symmetric  about  the  center  of  the  aperture  plane. 

This  results  in  the  sum  of  the  imaginary  part  of  the  P  component  equaling 
zero  under  all  undisturbed  conditions,  where  undisturbed  is  defined  as 
a  perfect  reflector  with  no  discs.  Note  also  that  the  cosine  term  in 
the  denominator  of  Eq  (14)  causes  the  zones  to  be  perpendicular  to  the 

line,  thus  the  zones  rotate  with  choice  of  far  f.^.J  i. 

The  real  part  of  the  P  component  therefore  controls  the  pattern 
entirely  in  the  undisturbed  case.  As  9  goes  to  zero,  k  goes  to  Infinity 
allowing  the  summing  of  the  unmodulated  electric  field  over  the  entire 
aperture  plane.  As  9  increases  (from  zero)  the  cosine  term  increasingly 
modulates  the  electric  field  until  the  sum  of  the  positive  and  negative 
sections  reaches  a  minimum,  producing  the  first  null.  As  0  further  increases, 
the  magnitude  of  the  sum  of  the  P  components  increases  until  the  maximum 
of  the  first  sidelobe.  Figure  3  is  approximately  the  maximum  of  the 
first  sidelobe  for  aperture  distribution  shown.  The  sum  continues  to 
oscillate  as  0  increases,  producing  the  characteristic  pattern.  The 
cosine  squared  term  in  Eq  (5)  has  been  neglected.  Its  contribution, 
however,  is  minimal  in  the  ranges  of  interest  (3  dB  at  45°). 

Using  this  simplified  model  the  effects  of  the  discs  can  be  examined. 

The  discs  in  the  model  produce  a  phase  shift  and  it  is  easily  seen  how 
this  may  cause  nulls.  Dominant  zone  is  defined  as  a  zone  in  which  the 
modulated  electric  field  is  the  same  sign  as  the  sum  over  the  entire 
aperture.  By  locating  the  disc  in  a  dominant  zone  such  that  it  produces 
a  180°  phase  shift  in  the  modulated  field,  the  area  affected  by  the  disc 


will  effectively  change  signs.  If  the  disc  is  the  right  size,  this  will 
cause  the  sum  of  the  modulated  electric  field  (the  real  part  of  the  P  ) 
to  equal  zero.  A  larger  disc  with  slightly  less  (or  more)  than  exactly 
180°  phase  shift  may  modulate  the  electric  field  to  produce  a  zero  sum. 

The  zeroing  of  the  real  part  of  the  P^,  however,  creates  a  problem 
since  the  symmetry  of  the  sine-modulated  electric  field  (the  imaginary 
part  of  P  )  has  been  destroyed — thus  creating  a  non-zero  sum.  This 
explains  why  it  requires  two  discs  to  produce  a  null  in  the  pattern.  In 
Jacavanco's  original  demonstration  of  this  procedure  he  used  two  discs 
on  the  same  side  of  the  aperture  (Jacavanco,  1984a).  However,  this  model 
shows  the  "optimum"  disc  placement  in  terms  of  using  the  smallest  discs 
for  minimal  disturbance  of  the  overall  pattern,  as  symmetric  with  respect 
to  the  center  of  the  dish  on  the  $'=$  line.  Using  this  placement,  the 
imaginary  part  of  the  P^'s  regain  their  symmetry  and  once  again  sum  to 
zero,  producing  a  zero  overall  sum  and  a  null  in  the  pattern. 

The  above  discussion  overlooks  one  problem  that  is  easily  solved 
but  requires  mention.  Even  the  simplified  disc  model  used  here  does  not 
produce  a  constant  phase  shift  over  the  area  affected.  (The  same  problem 
arises  using  over-sized  discs  adjusted  to  less  than  one-half  wavelength 
phase  shift.)  This  lack  of  an  exact  one-half  wavelength  phase  shift  is 
not  a  problem  in  zeroing  of  the  real  part  of  the  P^.  However,  even  with 
symmetrically  placed  discs  as  described  above,  a  non-symmetric  phase  shift 
again  causes  a  non-zero  imaginary  part  of  P^.  A  solution  is  easily  obtained 
by  slight  adjustment  of  the  phase  shift  produced  by  the  discs  to  produce 
a  zero  imaginary  and  real  part  of  the  P  component  (providing  the  discs 
are  large  enough) . 


Dominant  Zones 


The  sign  of  the  sum  of  the  modulated  electric  field  (real  part  of 
the  P  component)  determines  which  of  the  zones  (positive  or  negative) 
on  the  aperture  dominate.  Thus,  by  appropriate  phase  shifting  (and  thus 
changing  sign)  over  a  portion  of  the  dominant  zone,  nulling  may  be  accom¬ 
plished.  Unfortunately,  the  determination  of  whether  the  sum  for  a  given 
far  field  point  is  positive  or  negative  is  not  straightforward. 

If  there  is  no  blockage  for  the  feed,  the  sign  of  the  sum  alternates 
at  each  sidelobe.  The  area  of  zero  electric  field  caused  by  the  feed 
blockage,  however,  destroys  this  simple  property  and  requires  a  more 
rigorous  approach. 

The  feed  blockage  may  be  viewed  as  a  negative  electric  field  super¬ 
imposed  on  the  aperture  in  the  blocked  region,  thus  effectively  producing 
zero  electric  field.  This  results  in  a  far  field  pattern  that  is  the 
superposition  of  two  Bessel  functions  (Silver,  1949:191).  The  first 
function  is  positive  with  a  narrow  main  beam  associated  with  the  entire 
aperture:  the  second  is  negative  with  a  very  broad  main  beam  associated 
with  the  blockage.  Note  that  the  Bessel  function  is  the  analytic  counter 
to  the  sum,  therefore  where  the  Bessel  function,  or  sum  of  Bessel  functions, 
is  positive,  the  sum  is  positive,  and  where  it  is  negative,  the  sum  is 
negative.  In  the  main  lobe,  the  positive  Bessel  function  dominates, 
and  the  sum  is  positive.  However,  the  main  lobe  of  the  small  aperture 
(the  negative  pattern  produced  by  the  negative  electric  field,  or  blockage), 
while  small  in  comparison  to  the  main  lobe  of  the  large  aperture,  may 
not  be  small  in  comparison  to  the  sidelobes  of  the  large  aperture.  Thus, 
over  the  area  outside  the  main  lobe  of  the  large  aperture,  but  within 


the  main  lobe  of  the  small  aperture  (the  blockage),  the  sum  still 
oscillates  but  remains  negative  due  to  the  dominance  of  the  smaller 
aperture.  (Note  that  this  analysis  is  predicated  on  approximating  the 
aperture  distribution  as  parabolically  tapered — see  Silver.) 

Disc  Placement 

There  are  two  methods  of  determining  the  dominant  zones  of  a  given 
far  field  point.  The  aperture  distribution  may  be  approximated  and  the 
sign  of  the  sum  determined  analytically  as  just  described,  or,  a  computer 
analysis,  such  as  the  one  in  the  appendix,  may  be  run  and  the  sign  of 
the  sum  determined. 

Once  the  dominant  zones  are  determined  the  question  of  in  which 
zone,  and  where  in  the  zone,  to  place  the  discs  arise.  If  the  aperture 
distribution  is  not  tapered,  the  answer  is  apparent  by  looking  at 
Figures  3,  4  and  5.  A  disc  centered  in  any  dominant  zone  (centered  in 
the  horizontal  sense,  as  shown)  will  couple  the  same  amount  of  power 
and  thus  produce  the  same  result.  However,  when  a  tapered  aperture  distri¬ 
bution  is  considered  (as  almost  all  realistic  distributions  are)  another 
factor  must  be  considered.  Obviously,  the  closer  to  the  center  of  the 
aperture  the  disc  is  located  the  more  power  is  available.  The  taper, 
therefore,  makes  the  dominant  zone  closest  to  the  center  of  the  aperture 
the  optimum  for  disc  placement.  The  disc  should  be  centered  vertically 
(see  Figure  3)  in  the  zone  (again,  the  closest  available  position  to 
the  center  of  the  aperture) .  The  horizontal  position  of  the  disc  in 
the  zone  will  not  be  centered,  but  slightly  offset  toward  the  center 
of  the  aperture  due  to  the  taper.  Compare  Figures  4  and  5  with 
Figures  6  and  7.  These  figures  are  all  the  same  aperture  and  far  field 


point,  with  the  only  difference  being  Figures  4  and  5  represent  a  uniformly 
illuminated  aperture  and  Figures  6  and  7  represent  a  parabolically  tapered 
aperture  distribution  on  a  20  dB  pedestal.  No  easy  solution  for  optimum 
horizontal  placement  has  been  found.  A  good  approximation  is  to  center 
the  disc  at  the  maximum  of  the  modulated  electric  field  (taking  into 
account  both  the  tapered  electric  field  and  sinusoidal  modulation). 

Disc  Size 

While  disc  placement  is  somewhat  straightforward,  the  required  disc 
size  is  not.  The  required  size  depends  on  (1)  depth  of  null  required 
(amount  of  cancellation) ,  (2)  aperture  distribution  produced  by  the  feed, 
and  (3)  disc  placement. 

The  depth  of  the  null  or  cancellation  is  directly  proportional  to 
the  required  disc  size.  Any  size  disc  properly  placed  will  produce  some 
cancellation.  In  this  paper  the  depth  of  null  desired  will  be  assumed 
to  be  the  noise  floor.  For  plots  and  computer  simulations,  80  dB  below 
the  undisturbed  main  beam  gain  will  be  the  assumed  noise  floor.  Therefore, 
the  size  of  the  discs  will  be  proportional  to  the  undisturbed  pattern 
level.  Enough  power  must  be  coupled  by  the  discs  to  cancel  the  undisturbed 
pattern. 

The  aperture  distribution  will  have  a  strong  effect  on  the  required 
disc  size.  The  amount  of  energy  coupled  by  a  given  disc  is  controlled 
by  the  energy  density  incident  upon  it.  For  a  highly  tapered  distribution, 
the  energy  density  is  very  high  near  the  center  of  the  aperture  and  thus 
a  smaller  disc  near  the  center  of  the  aperture  can  couple  as  much  energy 
as  a  larger  disc  in  the  same  place  with  a  less  tapered  distribution. 


The  effect  of  disc  placement  on  required  size  may  be  broken  into 


two  parts.  The  first  follows  from  the  preceding  paragraph.  If  the 
aperture  illumination  is  tapered,  the  closer  to  the  center  of  the  aperture 
the  disc  is  placed  the  more  energy  it  will  couple  for  a  given  disc  size 
(of  course,  the  disc  must  be  in  a  dominant  zone  to  produce  any  cancellation). 
As  discussed  in  the  previous  section,  disc  placement  in  the  zone  also 
plays  a  part.  Due  to  the  sinusoidal  modulation,  more  apparent  power 
is  available  near  the  center  of  the  zone  than  at  the  edge. 

Performance  of  an  Arbitrary  Disc  System 

All  discussion  to  this  point  has  assumed  that  the  discs  used  were 
small  enough,  and  placed  correctly,  to  be  contained  in  a  single  zone. 

In  practice  this  may  not  be  the  case.  The  position  of  the  desired  null 

^  may  not  be  known  prior  to  antenna  design,  or  it  may  be  desired  to  move 

the  null  during  operation.  The  question  is  then  the  capability  of  a 
given  disc  system  (size  and  location)  to  produce  a  null  at  an  arbitrary 
point.  This  is  similar  to  the  question  discussed  in  the  previous  section. 

The  difference  is  that  the  disc  may  not  be  entirely  in  one  zone,  thus 
adding  another  complication. 

The  capability  to  produce  cancellation  is  governed  by  the 
(1)  illumination  or  energy  density  incident  on  the  disc,  and  (2)  where 
the  disc  is  located  in  the  zones  at  a  given  angle.  (If  the  disc  is  not 
in  a  dominant  zone  It  may  not  be  able  to  produce  any  cancellation — only 
gain.)  If  the  disc  spans  more  than  one  zone,  the  area  of  the  disc  capable 
of  cancellation  is  reduced,  as  the  area  in  a  non-dominant  zone  effectively 
cancels  area  in  the  dominant  zone.  This  is  not  one-to-one:  the  effective 

t.  A 

cancelled  area  will  depend  on  relative  illumination  and  phase  shift  produced 
by  the  disc. 
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An  analysis  of  an  arbitrary  disc  system  follows  a  two  step  process. 

The  first  involves  determing  the  discs'  position  in  the  zones  produced 
by  the  desired  location  of  the  null  in  the  far  field.  The  second  step 
is  calculation  of  the  effective  area  of  the  discs.  The  effective  area 
of  the  disc  system  (area  of  the  disc  system  that  is  not  cancelled  by 
an  equivalent  area  in  a  non-dominant  zone)  directly  controls  the  amount 
of  cancellation  available. 

For  an  aperture  with  a  uniform  illumination,  effective  area  may 
be  roughly  approximated  by  the  disc's  physical  area  in  a  dominant  zone(s) 
minus  that  area  in  a  non-dominant  zone(s).  This  approximation  involves 
estimating  the  effect  of  the  modulation  as  a  factor  of  one-half  and 
cancelling  this  with  the  factor  of  two  available  from  the  disc's  ability 
not  to  just  zero  an  area's  contribution  (to  the  sum),  but  to  produce 
a  negative  of  that  contribution.  If  the  undisturbed  pattern  is  x  dB 
below  maximum  pattern  gain,  an  effective  area  of  the  aperture  may  be 
considered  to  cause  the  pattern  level  by  the  relationship. 

A  2 

x  =  10  log10  (15) 

a 


where  A  is  the  effective  area,  causing  the  pattern  level,  and  A  is  the  total 
e  a 

area  of  the  aperture.  An  effective  disc  area  equal  to  the  effective 

area  producing  the  sidelobe  is  necessary  to  produce  a  null.  Therefore, 

the  minimum  total  effective  disc  area  (A,)  is 

d 


Unfortunately,  when  a  realistically  tapered  aperture  field  is  considered 
this  approximation  breaks  down  and  a  more  rigorous  approach  is  required 
(a  numerical  evaluation  will  be  necessary  for  most  distributions) . 


Effects  of  a  Non- Ideal  Antenna 

All  discussion  up  to  this  point  has  assumed  an  ideal  antenna  analyzed 
using  only  geometrical-optics  "modes."  While  this  approach  certainly 
simplifies  the  model  to  a  reasonable  level  of  complexity,  questions  arise 
as  to  the  validity  of  certain  assumptions.  Therefore,  a  short  discussion 
of  these  assumptions,  and  their  effects,  follows. 

Aperture  Integration  and  Geometrical  Optics .  Aperture  integration 
is  an  inexact  analysis  only  to  the  point  that  (1)  not  all  radiation  from 
the  antenna  passes  through  the  chosen  aperture,  and  (2)  the  actual  fields 
in  the  aperture  are  not  accurately  determined.  The  solution  for  the 
first  problem  is  to  choose  the  appropriate  aperture  such  that  the  far 
field  pattern  in  the  area  of  interest  is  accurately  determined.  For 
reflector  antennas  it  is  well-known  that  the  "ray-optics"  aperture  used 
in  this  model  gives  acceptable  results  within  45°-60°  of  the  main  beam 
(Silver,  1949:158). 

The  use  of  the  geometric  optics  field  in  the  aperture  is  obviously 
an  approximation.  While  this  field  is  present,  there  are  many  other 
contributors  to  the  total  field — scattering  from  feed  supports,  direct 
radiation  from  the  feed,  and  diffraction  from  the  edge  of  the  dish.  The 
assumption  made  here  is  that  these  fields  are  small  in  comparison  to 
the  geometric  optics  fields  and,  in  any  case,  are  subject  to  cancellation 
by  the  discs  in  the  same  manner  as  the  geometric  optics  field. 

Effects  of  Reflector  Imperfections .  One  of  the  major  problems  with 
reflectors  of  this  type  is  the  degree  to  which  the  reflector  dish  actually 
conforms  to  a  parabaloid  of  revolution.  This  model  represents  a  perfect 
surface;  however,  realistically  many  reflectors  of  this  type  may  have 


substantial  errors. 


The  effect  of  distortions  of  the  surface  may  be  modeled  in  Eqs  (12) 
and  (13)  as  an  additive  term  on  the  argument  of  the  sinusoids.  The  effect 
is  to  shift  the  zones  on  the  aperture.  If  the  distortion  is  very  small 
(in  terms  of  a  wavelength),  the  effect  may  be  disregarded.  However, 
if  the  distortion  is  a  significant  fraction  of  a  wavelength,  the  effect 
of  the  actual  dish  geometry  must  be  considered  in  any  analysis.  A 
distortion  of  one  quarter  wavelength  will  cause  a  phase  shift  of  180° 
and  thus  a  total  reversal  of  the  zone. 

Idealized  Feed.  The  effects  in  an  actual  antenna  system  of  having 
a  non-ideal  feed  may  be  separated  into  two  parts.  The  first  is  extraneous 
radiation  from  the  feed.  This  is  considered  in  the  previous  section 
concerning  the  additional  fields  in  the  aperture  not  due  to  geometrical 
optics. 

The  second  point  is  that  a  feed  may  not  produce  an  aperture 
distribution  that  is  parabolically  tapered.  A  review  of  the  analyses 
will  reveal  that  although  this  distribution  is  selected  as  mathemati¬ 
cally  convenient,  in  terms  of  nulling,  the  distribution's  effect  is 
limited  to  power  available  at  a  given  point  on  the  aperture.  Thus,  while 
the  actual  aperture  taper  may  affect  required  disc  size,  and  optimum 
disc  placement,  it  has  no  effect  on  the  basic  mechanism  of  nulling. 


Computer  Analysis 

The  equations  of  the  preceding  chapter  were  used  to  write  Fortran 
code  to  evaluate  the  effect  of  the  discs  on  the  far  field  pattern  of 
a  reflector.  This  code  was  used  to  test  the  model's  ability  to  actually 
produce  nulls  in  the  far  field  pattern  and  investigate  properties  of 
the  system.  The  code  and  a  brief  explanation  of  its  use  is  in  the  appendix; 
no  attempt  will  be  made  to  explain  its  working  here. 

System  Modeled 

The  antenna  modeled  is  a  four-foot  diameter  parabolic  reflector 
with  a  focal  length  of  27  inches.  It  operates  at  2.8  GHz  and  is  fed 
to  produce  a  parabolic  aperture  distribution  (order  one)  on  a  20  dB  pedestal. 
The  feed  polarization  is  assumed  to  be  purely  linear  (y) .  This  produces 
an  aperture  distribution  of  the  same  polarization  except  as  corrupted 
by  the  dish  (Stutzman  and  Thiele,  1981:427).  The  feed  blockage  is  assumed 
to  be  a  six- inch  diameter  circle  at  the  center  of  the  aperture.  These 
parameters  were  chosen  to  model  a  generic  antenna. 

Two  disc  systems  are  modeled.  The  first  (disc  system  I)  has  two 
seven-inch  diameter  discs  symmetrically  centered  9.5  inches  from  the 
center  of  the  aperture  on  the  x-axis  of  the  aperture  plane.  The  second 
(disc  system  II)  has  2. 5- inch  discs  similarly  placed,  but  centered  4.5 
inches  from  the  center  of  the  aperture  plane. 

These  two  disc  systems  were  selected  because  they  demonstrate  not 
only  the  systems'  ability  to  produce  nulls,  but  also  the  range  and  limita¬ 
tions  of  this  system.  Note  also  that  the  discs  do  not  overlap  and  thus 
could  both  be  mounted  on  the  dish  simultaneously  to  enhance  system  flexibility 
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Radiation  Patterns 


All  the  radiation  patterns  shown  are  H-plane  patterns.  The  reason 
for  this  choice  is  that  both  disc  systems  are  optimally  placed  to  produce 
nulls  in  the  <p-0  (H)  plane.  (The  preceding  chapter  explains  the  reasons 
a  disc  system  symmetrically  placed  across  the  center  of  the  aperture 
on  the  line  is  optimum.) 

The  first  pattern  shown  (Figure  8)  is  the  radiation  pattern  produced 
by  the  antenna  with  no  discs  mounted.  This  is  shown  for  two  reasons. 

The  first  is  to  demonstrate  that  the  code  used  will  accurately  implement 
the  equations  of  chapter  one.  The  radiation  pattern  of  the  antenna  without 
the  discs  may  be  solved  analytically,  and  the  agreement  between  the  analytic 
solution  and  the  numerical  solution  shown  is  excellent  (within  . 1  dB) . 

The  second  reason  for  Figure  8  is  to  provide  a  baseline  for  comparison. 
Important  features  to  note  from  this  pattern  are:  a  6°  beamwidth,  an 
approximately  -21  dB  first  sidelobe,  and  sidelobe  maxima  at  10°,  20° 
and  30°. 

Disc  System  _T 

The  parameters  of  this  system  were  selected  with  two  goals.  First, 
a  null  is  to  be  placed  at  10° — the  maximum  of  the  first  sidelobe.  The 
first  sidelobe  is  the  most  difficult  position  to  null  because  it  has 
the  most  power.  This  may  also  make  it  the  most  important  position  to 
be  able  to  null  in  an  actual  system.  The  second  consideration  was  to 
make  the  system  as  flexible  as  possible  (flexible  in  terms  of  being  able 
to  null  over  a  large  range  of  9). 

The  second  specification  makes  the  optimal  position  of  the  preceding 
chapter  no  longer  the  best.  By  mounting  discs  closer  to  the  center  of 
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the  aperture  than  the  optimal  position,  nulling  is  possible  over  a  larger 
range  of  0.  Our  main  interest  is  the  ability  to  produce  nulls  as  0  grows 
larger;  no  real  attempt  is  made  to  null  into  the  main  beam.  The  cost 
of  this  flexibility,  of  course,  is  that  a  larger  disc  is  required  in 
this  position  than  would  be  in  the  optimal  position.  The  performance 
of  this  system  is  shown  in  Figures  9  through  16. 

The  fulfillment  of  the  first  requirement  is  shown  in  Figure  9.  By 
proper  adjustment  of  the  two  discs,  a  null  is  produced  at  10°. 

Figures  10  and  11  show  the  effects  of  the  discs  on  the  aperture 
electric  field  along  the  x-axis  (discs  extend  from  -13  and  -6,  and  6  to 
13  inches).  Note  the  approximate  180°  phase  shift  over  the  affected  area 
characteristic  of  a  minimum  sized  disc. 

This  disc  system  is  capable  of  producing  a  null  in  the  radiation 
pattern  anywhere  from  7°  to  20°.  Figure  12  shows  a  null  placed  at  20° 

(the  second  sidelobe  of  the  undisturbed  pattern).  Figures  13  and  14 
show  the  discs'  effect  on  the  aperture  electric  field  is  no  longer 
contained  in  one  zone.  (For  all  sidelobes  out  to  45°  the  negative  zones 
are  dominant  for  this  antenna  system.)  Also,  the  phase  shift  produced 
by  the  discs  is  much  less  than  180°. 

The  extremes  of  the  capabilities  of  this  disc  system  are  shown  in 
Figures  12  and  16.  The  null  shown  in  Figure  16  at  7°  is  the  "inner  limit" 
of  the  nulling  capability  (within  a  whole  degree).  Note  that  this  is 
actually  inside  the  first  null  of  the  undisturbed  pattern.  Figure  12 
shows  the  system  adjusted  to  produce  a  null  at  20°,  the  "outer  limit" 
of  this  system's  capability. 
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Disc  System  II 


The  parameters  for  this  disc  system  were  selected  to  allow  the  use 
of  the  smallest  discs  possible  to  produce  a  null  which  was  at  31°  (The 
maximum  of  the  third  sidelobe) .  It  was  desired  to  have  the  effective 
nulling  range  of  this  system  not  overlap  with  that  of  system  I  in  order 
to  demonstrate  the  ability  of  a  set  of  disc  systems  to  null  where 
neither  is  totally  effective  individually. 

Disc  system  II,  adjusted  to  produce  a  null  at  31°,  is  shown  in 
Figure  17.  The  modulated  electric  fields  in  this  configuration  are 
shown  in  Figures  18  and  19.  Again,  note  the  approximate  180°  phase 
shift  over  the  area  affected  by  the  disc. 

This  system  is  capable  of  generating  nulls  from  24°  to  in  excess 
of  60°  which  is  the  outer  range  where  this  model  is  considered  valid. 

The  radiation  pattern  of  the  antenna  with  the  disc  system  adjusted  to 
produce  a  null  at  24°  is  shown  in  Figure  20. 

A  null  in  the  radiation  pattern  at  45°  is  shown  in  Figure  21.  This 
pattern  is  of  interest  for  two  reasons.  Examining  Figures  22  and  23, 
the  modulated  electric  fields  along  the  x-axis  for  this  configuration 
reveal  the  narrow  width  of  the  zones  at  this  angle  (one  quarter  wavelength) . 

A  second  interesting  property  of  this  nulling  method  becomes  evident 
comparing  Figures  17  and  21.  In  Figure  17  the  desired  null  is  at  31° 
and  the  pattern  at  45°  is  approximately  at  the  same  level  as  the  undisturbed 
pattern.  Figure  21  shows  the  desired  null  at  45°,  but  also  shows  the 
pattern  level  at  31°  is  over  40  dB  down  from  the  undisturbed  pattern. 

Thus,  Figures  17  and  21  demonstrate  that  one  disc  system  may  be  able 
to  provide  a  reasonable  amount  of  cancellation  to  more  than  one  noise 
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source  if  both  sources  were  considered  at  the  same  time.  (This  was  not 
feasible  using  the  simple  code  in  the  appendix;  however,  in  an  actual 
system  the  algorithm  would  likely  have  to  respond  to  total  noise  power 
because  no  distinction  would  be  available  between  noise  from  different 
angles . ) 

Disc  System  III 

What  is  now  defined  as  disc  system  III  is  actually  the  combination 
of  systems  I  and  II.  As  noted  earlier,  these  two  systems  do  not  overlap, 
and  thus  could  be  mounted  on  the  same  dish.  The  problem  we  wish  to  address 
is  whether  the  two  systems  together  can  accomplish  what  neither  can  indepen¬ 
dently,  i.e.,  provide  a  nulling  capability  between  20°  and  24°.  Figure  25, 
with  a  null  at  22°,  demonstrates  that  this  is  indeed  possible;  the  effects 
of  the  two  systems  may  be  seen  to  add  directly.  This  would  be  expected 
since  the  discs  simply  act  to  reduce  the  magnitude  of  the  P-component. 

This  effect  adds  in  a  linear  fashion;  however,  since  the  pattern  strength 
is  proportional  to  the  magnitude  of  the  P-coraponent  squared,  the  amount 
of  cancellation  available  is  much  greater  than  a  linear  sum.  This  model 
shows  both  disc  systems  must  be  able  to  provide  some  cancellation;  if 
a  system  is  not  able  to  provide  any  cancellation,  adding  it  to  another 
system  simply  degrades  the  performance  of  the  effective  system. 
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Figure  17.  H-Plane  Radiation  Pattern  With  Null  at  9-31°  (Disc  System  II) 
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Figure  23.  Sine  Modulation  on  x-Axis  With  Null  at  9=45°  (Disc  System 
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Figure  25.  H-Plane  Radiation  Pattern  With  Null  at  9=22°  (Disc  System  III) 
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System  Bandwidth 

The  other  effect  investigated  using  this  model  was  the  effect  of 
frequency  on  the  antenna  pattern  containing  nulls.  Since  most  noise 
sources  contain  more  than  a  single  frequency  (except  perhaps  a  tone  jammer) , 
the  width  of  a  null  in  frequency  has  a  definite  bearing  on  system  capability 
in  any  realistic  situation.  The  feed  (or  receiver)  in  an  actual  system 
may  provide  noise  rejection  outside  the  design  bandwidth.  However,  most 
systems  must  be  designed  to  operate  over  some  finite  bandwidth  and  are 
susceptible  to  noise  in  this  bandwidth. 

Frequency  Performance  Plots 

The  disc  systems  used  to  generate  these  plots  are  the  same  as  used 
in  the  previous  sections.  Referring  to  Figures  26  through  29,  actual 
pattern  points  were  calculated  only  at  the  frequencies  shown  as  asterisks 
on  the  "without  discs  plot."  Thus,  the  data  is  actually  only  exact  at 
the  thirteen  frequencies  with  a  straight  line  connecting  these  points 
to  show  an  estimate  of  performance. 

A  point  to  consider  analyzing  these  plots  is  that  while  performance 
is  examined  over  a  20%  bandwidth,  the  discs  were  adjusted  to  provide 
cancellation  at  the  center  frequency.  The  code  used  has  no  provision 
for  cancellation  over  a  finite  bandwidth.  An  actual  system,  however, 
could  be  designed  to  do  so.  While  no  claim  is  made  that  the  discs  could 
be  adjusted  to  provide  a  -80  dB  null  over  a  finite  bandwidth,  a  better 
adjustment  than  shown  may  be  possible  (in  terms  of  minimizing  total  noise 
power  to  the  receiver)  if  the  adjustment  were  done  with  respect  to  the 


entire  bandwidth. 


Figure  26.  This  plot  shows  the  performance  at  7°  with  disc  system  I 
adjusted  in  the  same  manner  as  Figure  16.  This  demonstrates  the  problem 
with  locating  a  null  in  close  proximity  to  a  null  in  the  undisturbed 
pattern  (pattern  of  the  antenna  without  discs).  The  change  in  frequency 
has  effectively  moved  the  first  null  in  Figure  8  "in"  to  7°  at  3.0  GHz. 
Thus  the  disc  system  provides  cancellation  performance  only  at  frequencies 
less  than  approximately  2.9  GHz. 

Figure  27.  This  is  the  performance  of  the  same  antenna  and  disc 
as  shown  in  Figure  10  at  10° .  While  performance  is  obviously  best  at 
the  center  frequency,  there  is  cancellation  over  approximately  a  10% 
bandwidth.  As  discussed  before,  it  may  be  possible,  with  different  disc 
adjustment,  to  provide  better  cancellation  over  the  entire  bandwidth 
by  sacraficing  some  cancellation  at  the  center  frequency. 

Figure  28.  This  is  disc  system  II' s  performance  at  20°  with  the 
same  adjustment  as  Figure  13.  This  shows  good  performance  over  almost 
a  20%  bandwidth  with  especially  good  performance  over  the  lower  half 
of  the  frequency  range. 

Figure  29.  This  figure  shows  the  performance  of  disc  system  II 
with  the  same  disc  settings  as  Figure  18  at  31°.  Again,  excellent  perfor¬ 
mance  is  noted  over  a  10%  bandwidth  with  cancellation  available  over 
almost  all  of  the  20%  bandwidth. 
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Conclusion 


The  theoretical  model  presented,  while  not  producing  any  exact 
equations  for  the  disc  system,  does  offer  several  "rules  of  thumb."  Any 
exact  formulation  would  have  to  be  suspect  regardless,  since  this  is 
only  a  first  order  model.  Even  assuming  an  ideal  feed  and  perfect 
reflector,  second  order  effects,  such  as  diffraction,  become  significant 
in  the  sidelobe  region. 

Disc  Effects 

The  discs  were  found  capable  of  generating  the  desired  nulls  in 
the  radiation  pattern  while  leaving  the  main  beam  gain  relatively  unaffected. 
The  maximum  gain  loss  is  only  slightly  over  one  dB  (1.2  dB  on  the  pattern 
with  a  null  at  10°). 

The  three  dB  beamwidth  was  affected  even  less,  although  as  pattern 
levels  were  calculated  only  every  degree,  exact  numbers  are  not  available. 

The  discs  do  cause  sidelobe  levels  to  be  raised;  however,  no  sidelobe 
was  raised  above  -20  dB  (only  1  dB  higher  than  the  undisturbed  first 
sidelobe  level).  Overall  performance  is  excellent,  and  in  general  agree¬ 
ment  with  performance  as  demonstrated  experimentally  by  Jacavanco 
(Jacavanco,  1984a). 

Dish  Size 

It  should  be  noted  that  all  work  done  (both  experimentally  by  Jacavanco 
and  in  this  model)  has  been  with  reflectors  that  are  relatively  small 
in  terms  of  wavelengths.  A  problem  is  that  while  the  zone  widths  are 
independent  of  the  reflector  size  (see  Eq  (14)),  required  disc  size  is 
not  (see  Eq  (16)).  Simply  enlarging  the  disc  in  direct  proportion  to 


the  reflector  therefore  will  not  give  equivalent  performance  with  a  larger 
reflector.  Because  the  larger  disc  will  span  more  zones,  its  effective 
area  will  be  smaller.  One  solution  would  be  to  use  an  increased  number 
of  smaller  discs  on  the  larger  reflector  (appropriately  placed  in  the 
zones).  Another  option  would  be  to  use  a  shape  other  than  the  circular 
disc  to  supply  an  equivalent  area  in  the  zone. 

Comparison  of  Model  with  Experimental  Results 

In  order  to  try  and  validate  this  model,  experimental  data  was  obtained 
from  Daniel  Jacavanco  (Jacavanco,  1984b).  All  attempts  to  recreate  or 
match  these  results  failed.  This  failure  is  thought  to  be  a  result  of 
poor  data  on  the  antenna  used  in  the  experiment,  not  a  failure  of  the 
basic  theory. 

Two  basic  characteristics  of  the  antenna  used  In  the  experiment 
were  not  well-known  when  the  attempt  to  match  results  was  made.  The  first 
is  the  actual  aperture  distribution  produced  by  the  feed.  The  distribution 
was  modeled  as  a  parabolic  taper  of  order  one  on  a  20  dB  pedestal.  Data 
obtained  on  the  radiation  pattern  of  the  feed  subsequently  has  shown 
that  although  the  pedestal  height  was  approximately  correct,  the  taper 
estimate  was  too  steep. 

The  second  error  in  the  antenna  as  modeled  was  the  assumption  that 
the  reflector  was  a  perfect  paraboloid.  The  actual  antenna  used  in  the 
experiment  was  a  general  purpose,  medium  performance  reflector.  When 
the  surface  of  the  reflector  was  measured  it  was  found  to  have  significant 
errors.  This  is  characteristic  of  medium  performance  reflectors,  but 


is  difficult  to  model. 


Either  of  these  errors  in  the  model  could  cause  failure  of  the  model 


to  match  experimental  results.  It  is  not  surprising,  then,  that  no  correla¬ 
tion  could  be  found  between  the  experimental  and  model  results.  Another 
attempt  is  being  made  to  correct  the  model  and  match  results.  This  should 
be  the  primary  goal  of  any  further  research. 


Areas  of  Further  Investigation 

There  are  several  areas  where  more  investigation  on  this  topic  may 
provide  valuable  results.  As  mentioned  earlier,  however,  any  further 
research  should  probably  begin  with  an  attempt  to  develop  a  more  accurate 
model  of  the  antenna  used  in  the  work  by  Jacavanco,  in  order  to  match 
experimental  results.  This  is  not  a  trivial  goal  as  errors  in  the  reflector 
surface  will  be  difficult  to  accurately  model.  To  establish  confidence 
in  this  approach,  however,  corroboration  by  experimental  results  is 
essential. 


Rectangular  Sections 
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The  zones  produced  by  this  approach  suggest  another  method  of  imple¬ 
menting  this  basic  idea.  Rectangular  sections,  instead  of  discs,  mounted 
with  their  long  axis  perpendicular  to  the  axis  of  the  desired  null  would 
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appear  to  be  a  more  efficient  system.  This  system  could  be  adjusted 
either  by  varying  the  sections'  distance  off  the  reflector  in  the  same 
manner  as  the  discs,  or  by  mounting  them  one-fourth  wavelength  off  the 


reflector  (to  produce  a  180°  phase  shift)  and  moving  them  laterally  to 
adjust  the  pattern.  Either  system  would  be  a  more  complicated  mechanical 
arrangement  than  the  discs,  but  performance  gains  may  be  worth  the  com¬ 


plexity. 


Polarization  dependence  is  a  problem  that  an  analysis  of  this  system 


would  have  to  consider.  If  the  sections  were  mounted  such  that  their 
long  axis  was  not  collinear  with  the  feed  polarization,  their  effective 
area  would  likely  be  reduced.  An  approach  other  than  geometrical  optics 
would  have  to  be  employed,  and  this  limitation  considered,  in  any  comparison 
of  performance. 

Adaptive  Array 

Another  area  of  interest  would  be  the  design  of  an  array  of  discs 
mounted  on  a  dish  to  provide  nulling  at  an  arbitrary  angle.  An  algorithm 
to  control  the  discs  to  provide  an  adaptive  system  could  then  be  considered. 
Main  beam  gain  would  suffer  excessively  if  too  many  discs  were  mounted 
on  the  dish.  However,  a  system  that  used  discs  made  of  metalized  plastic 
that  could  be  pulled  back  flush  to  the  dish  when  not  in  use  might  overcome 
this  problem.  An  algorithm  that  could  efficiently  control  such  an  array 
would  be  a  challenging  project.  If  such  a  system  could  be  shown  to  be 
feasible,  it  could  provide  an  attractive  alternative  to  a  conventional 
adaptive  array  in  some  applications. 

Number  of  Discs  Required  Per  Null 

The  array  problem  brings  to  light  another  question  that,  while 
theoretical,  is  of  practical  interest.  The  model  as  presented  in  this 
paper  demonstrates  that  the  generation  of  a  single  null  requires  two 
discs.  Daniel  Jacavanco  has  proposed  that  in  general  to  generate  n  nulls 
in  the  far  field  pattern  would  require,  following  from  array  theory, 
n+1  discs  (Jacavanco,  1984a).  The  theory  of  this  paper,  however,  would 
seem  to  suggest  that  two  discs  per  null  may  be  required,  as  for  each 


null  you  must  control  two  variables  (the  imaginary  and  real  parts  of 
the  P-component) .  Because  these  variables  are  not  independent,  however, 
it  is  not  immediately  clear  whether  two  discs  are  required  or  not.  An 
investigation  of  the  number  of  discs  required  to  produce  n  nulls  at  arbitrary 
locations  would  definately  be  of  interest  to  any  future  research  on  arrays 
of  discs. 

Optimum  Distortion 

The  discs  are  a  method  of  reflector  distortion  that,  while  easily 
fabricated  and  adjusted,  do  not  offer  an  optimum  solution.  Eqs  (10) 
and  (11)  do,  however,  offer  a  solution  for  an  optimal  distortion.  If 
the  distortion  induced  was  such  that  the  argument  of  the  sinusoids  was 
equal  to  an  even  or  odd  integer  multiple  of  pi  (depending  upon  which 
zone  was  dominant),  the  maximum  cancellation  per  area  distorted  would 
be  realized.  Note  that  two  mirror  image  areas,  symmetrically  placed 
across  the  aperture,  would  still  be  required  to  insure  a  zero  contribution 
from  the  imaginary  part  of  the  P-component.  Note  also  a  perfect  surface 
over  the  rest  of  the  reflector  is  required  to  insure  optimal  null 
generation. 

A  scheme  of  this  sort,  while  optimal  (in  terms  of  minimum  area 
distorted  per  cancellation),  may  not  be  preferred  in  an  actual  implemen¬ 
tation.  One  of  the  main  problems  with  this  solution  is  that  it  would 
require  the  exact  desired  null  position  to  be  known.  While  it  is  relatively 
easy  to  imagine  an  array  of  discs  controlled  by  some  algorithm  to  provide 
a  degree  of  adaptive  nulling,  no  simple  method  is  readily  apparent  to 
adaptively  implement  this  scheme. 


Another  problem  that  this  system  would  appear  to  have  is  that  since 
the  solution  is  optimal,  it  would  be  much  more  sensitive  to  variations 
in  either  angles  or  frequency.  Sensitivity  to  the  noise  sources'  movement 
in  the  far  field  may  or  may  not  be  a  problem  depending  on  scenario.  However, 
a  high  sensitivity  to  frequency  (i.e.,  a  narrow  bandwidth  null)  would 
cause  problems  in  almost  any  imaginable  implementation. 

Summary 

The  initial  objective  was  to  develop  a  model  that  would  explain 
the  nulling  performance  of  discs  mounted  on  a  reflector  as  demonstrated 
by  Jacavanco  (Jacavanco,  1984a).  Aperture  integration  was  selected  as 
a  well-known  technique  that  appears  well-suited  to  the  analysis  of  this 
antenna  system.  While  the  model  developed  is  only  a  first  order  approxi¬ 
mation,  the  important  aspects  of  the  system  appear  to  be  accurately  repre¬ 
sented. 

The  model  developed  is  broken  down  mathematically  to  provide  insight 
into  the  mechanism  of  null  generation  and  its  control.  Guidelines  are 
developed  for  disc  placement  for  null  generation  based  on  this  model. 

No  satisfactory  guidelines  were  found  for  required  disc  size  although 
a  rough  method  of  estimation  is  presented. 

The  mathematical  model  is  implemented  in  Fortran  code  and  is  exercised 
to  demonstrate  its  ability  to  generate  nulls.  The  computer  model  exhibits 
the  properties  predicted  by  the  basic  theory  and  is  used  to  further  examine 
system  characteristics  such  as  bandwidth  performance.  Theoretical  results 
from  the  model  were  compared  to  experimental  results  without  success. 

This  is  thought  to  be  a  result  of  too  idealistic  an  antenna  model  rather 
than  error  in  the  basic  theory. 
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Appendix 

Computer  Codes 

Listed  are  the  codes  used  in  the  computer  analysis  of  this  work. 

The  codes  are  a  straightforward  implementation  of  the  equations  of 
Chapter  One.  Both  are  written  in  standard  Fortran  77. 

Method 

The  aperture  is  partitioned  into  hexagons,  and  the  hexagon  sampled 
at  7  points  to  approximate  the  double  integral  (see  Abramowitz  and  Stegun, 
1972:894).  The  variables  PXR,  PXI,  PYR,  and  PYI  are  the  variables  of 
summation  for  the  real  and  imaginary  part  of  the  x-  and  y-component  of 
the  P  vector.  These  variables  are  then  used  to  calculate  pattern  strength 
that  is  normalized  to  the  pattern  strength  at  0=  4>=0 . 

Assumptions 

There  are  several  assumptions  built  into  the  code  which  must  be 
recognized.  The  first  is  that  this  code  models  a  parabolic  reflector. 

This  is  basic  to  the  program  and  is  not  easily  changed. 

Several  assumptions  are  made  concerning  the  feed.  It  is  assumed 
to  block  a  circular  area  (variable)  at  the  center  of  the  aperture  and 
produce  a  purely  y-polarized  field  (the  coordinate  system  used  in  the 
computer  code  is  the  same  as  in  the  rest  of  this  paper) .  Neither  of 
these  parameters  may  be  changed  by  specification  from  the  files.  The 
feed  is  also  assumed  to  produce  a  parabolically  tapered  electric  field. 

A  wide  range  of  aperture  distributions  may  be  approximated  using  this 
formula  and  it  is  easily  changed.  Any  function  desired  may  be  substituted 


in  the  EFIELD  function. 


The  code  accesses  the  files  as  listed  following  the  two  programs 
to  describe  the  antenna  modeled  and  the  far  field  points  of  interest. 

The  file  ANTENNA  lists  all  the  parameters  of  the  basic  antenna,  and  the 
file  DISCS  provides  the  parameters  of  the  disc  system.  ANGLES  controls 
the  far  field  points  of  interest  and  the  number  of  these  points. 

The  two  programs  differ  in  that  MA1N1  reads  the  angles  theta  from 
the  file  ANGLES  and  outputs  both  theta  and  the  values  of  the  summation 
variables  to  the  file  PS.  Program  MAIN2  reads  the  values  of  theta  from 
the  file  PS  along  with  the  variables  of  summation.  When  program  MAIN1 
is  run  first  for  a  given  antenna  system  (and  angle  phi),  it  loads  the 
file  PA  with  the  appropriate  values.  Program  MAIN2  may  then  be  run  with 
different  disc  configurations.  Since  MAIN1  integrates  over  the  entire 
aperture  and  MAIN2  only  over  the  discs  (using  the  values  of  the  summation 
variables  to  calculate  the  value  of  the  integral  over  the  rest  of  the 
aperture),  MAIN2  requires  a  much  shorter  run  time.  The  code,  as  written, 
is  a  brute  force  integration  technique,  with  no  thinning  as  is  commonly 
employed  in  this  type  of  program.  (This  was  done  to  insure  the  discs' 
effects  were  properly  calculated.)  Program  MAIN2  is  an  attempt  to  provide 
a  much  more  efficient  analysis  tool.  For  the  antenna  analyzed  in  this 
paper  (the  parameters  of  which  are  shown  in  the  listing  of  the  files 
in  this  appendix),  MAIN1  required  approximately  8  cpu  seconds  to  run 
on  a  CDC  Cyber  750  (a  similar  6-foot  dish  required  approximately  16  cpu 
seconds).  MAIN2,  however,  required  only  slightly  over  1  cpu  second  to 
run  on  the  same  machine  with  the  two  seven-inch  discs  (note  that  run 
time  for  MAIN2  is  controlled  by  disc  size;  a  6-foot  dish  would  require 


the  same  amount  of  time  as  a  4-foot  dish  if  the  area  of  the  disc  systems 
of  the  two  were  equal) . 

Both  programs  output  the  pattern  levels  to  the  file  PATTERN,  and 
to  the  terminal.  (This  assumes  the  programs  are  being  run  interactively 
if  this  is  not  the  case,  the  deletion  of  the  one  output  line  will  cancel 
the  terminal  output.) 


PROGRAM  MAIN1 

D I  MENS  ION  W t7) *DX  t7 ) *DY  <7>  .UAEF  <2>*  XDt  4>  , 

*YOt 4) *8Pt 4>  »OSHIFTt 4>  *  DISCS (4) 

SEAL  MNt4) 

D0U8LE  PRECISION  PXR *PX I *PYR *PY I 
0PENt7,FILE  =  « ANTENNA*  > 

0PENt8*FILE=*  ANGLES  ' ) 

0PENt9,FILE=*DISCS*  ) 

DPENtl3,FIL£=*PATTERN«> 

OPENtll *FILE=  »PS»  ) 

DATA  DX ( 1 ) *  DY 1 1 ) »  DY (2 ) • OYt 3 ) /4* 0 . C/ 

DATA  PXR»PXI*PYR*PY 1/4*0. 0/ 

REWIND! 7 ) 

REU I NDt  8 ) 

REUINO 1 9  > 

C  ANTENNA  PARAMETERS  READ 

READ(7»10  0C )SI D£« RADIUS  *F RE Q*Ft  BLOCK *C  »N,PMAX  *IDISC 
10  00  FORMATt/* 6t/*E14.7* /) f /1 2 « / « /« E 1 4 .7 */ « / « 1 3 ) 

C  FAR  FIELD  PLANE  OF  INTEREST  CPHI)  AND  NUMBER  OF 

C  POINTS  DEFINED 

READt8.2Q0C  >PHI ,NANG 
2000  FORMAT  t/»/»F10.7*/»/*I3#/) 

C  DISC  PARAMETERS  READ*  AND  CONSTANTS  CALCULATED 

DO  50  IR=1* IDISC 

RE AD <9, 2600 )XD<IR)*YO(IRI *D ISCR (I R ) .DSHIFTtlR) 

2600  F0RMATt3t/*£14.7,/> ,/*£14.7> 

OELT A=DSHIFT tIRl *3. GE* 08/FREQ 
RP=SQRTtXD<IR)**2*YDt IRI**2> 

MNt IR)=-RP/t2.0*F  ) 

8PtIR)=tRP**2*4.0*F**2)/(2.0*RP)- 
♦DELTA*t2.0*F-MNtIR> *RP ) /( RP *SQR T ( 1. C* MNt IR I ** 2 ) ) 

50  CONTINUE 

C  CONSTANTS  FOR  INTEGRATION  SCHEME  ARE  ASSIGNED 

D I E= 3. 141 59 26535898 
‘3ETA=2.0*PIE*FREQ/3.0E+08 
PHI=PHr*PIE/180  .0 
Wtl>=258.  0/1008.0 
Wt2  >  =  125.0/1008.0 
U  t  3 ) =W 1 2 ) 

Wt4)=Wt2) 

Mt 5  >=Wt2> 

Wt6)=WC2> 

Ut  7>  =  Wt2> 

OX t2)=SIDE*SQRT tl4. 0>/5.0 
DX  t  3 >  =  -DX  1 2 > 

DXt5)=SIDE*SQRTtl4.  O/IO.O 
DX  t  7)=DX  t  5) 

DX  t  4>  =  -DX  t5 > 

DX  t  6 ) =DX  t  4 ) 

DYt  4)=SIDE*SQRTt42.0»/10.0 


o  y  <  5) =oy< a> 

or  {it)  =  -oy(4) 

DY ( 7>=OY<  6) 

C PH I ~ COS  <  PH  I) 

SPHI  =  SIN<PHI) 

Y INCR=3.0*SIDE 
XINCR=SaRT<3.0)+SIDE 
N1=INT (RAOIUS/XINCR) 

FAR  FIELD  POINTS  DEFINED  CTHETA) 

DO  700  L=ltNANG 
READ(8  *  3000  )THETAD 
FORMAT  (F7 .3 ) 

THETA=THETAD*PI£/180.0 
CTHETA=COS  CTHET  A) 

STHETA=SIN<  THET  A) 

XC=XINCR/2.0 

SUMMATION  OVER  THE  APERTURE 
DO  600  IX=1*2*N1*1 
XC=-XC*XINCR/2.0 
IF  CIX.EQ.1)XC=0.0 
I  T  E S  T  =  MOO ( I X  *  2  > 

IF  CITEST.GT.C)  THEM 

YC=0.0 

YCC=0.0 

GO  TO  300 

ELSE 

YC=YINCR/2.0 
YCC=YC 
GO  TO  300 
EMDIF 

YC=YC+YINCR 
GO  TO  300 
XC=-XC 

yc=ycc 

GO  TO  300 
YC  =  -YOYINCR 
CONTINUE 
IDTEST=0 

°Q  I  NT  IS  TESTED  TO  DETERMINE  IF  IT  IS  AF 
00  350  IO=l»IDISC 

DIST=SQRT<  (XC-XD(ID) )**2*<YC-YD(ID))**2 
IFCDIST.LT. DISCR(IO))  GO  TO  375 
CONTINUE 
I D  TEST =2 
CONTINUE 

RPRIME=SQRT(XC*XC+YC*YC> 

POINT  IS  TESTED  TO  DETERMINE  IF  IT  IS  ON 
IF  CRPR IME.LT. BLOCK )  GO  TO  100 
IF  (RPRIME.GT .RADIUS)  GO  TO  550 


AFFECTED 


DISC 


APERTURE 


c 


I  NO  IV  I  DUAL  ELEMENT  13  SUMMED  OVER 
00  500  1=1*7 
<  =  XOOV(I  ) 

Y  =YC ♦OY  ( I > 

RPR  I  ME=SQ.RT  (  X  *X  +Y  *  Y  ) 

3HIP=ATAN2( Y,X> 

THETAP=2.0»ATAN(RPRIME/(2.0*F) I 
CPP=CDS (PHIP) 

SPP=SI"J(PHIP) 

CTP=COS(TH£TAP> 

STP=SIN(THETAP> 

O=SQRT(1.0-SPP*SPP*STP*STP) 

UAEF(1)=( SPP*CPP* ( 1 • 0“C  TP7/D) 

UAEF(2)=( SPP*SPP*CTP+CPP*CPP)/D 
C2=EFIELD (RPR IME  * RADIUS  »C,N  >*W(I > 

ARG=8ETA*RPRI ME*STH  ET  A*COS( PHI-PH IP ) 

IF  (IOTEST.EQ.O  )  THEN 
XP=(RPRIME-8P(ID> )*  MN(ID) 

3  =  SQRT(XP**2«-RPRIME**2)«-XP 
SMIFT=(2.0*F-S)*BETA 
ARGO=ARG- SHIFT 
PR=C2*C0S  (ARGO) 

3I=C2*SIN(ARGD) 

ELSE  ‘ 

PR=C2*C0S (ARGI 
3I=C2*SIN(ARG) 

ENOIF 

PXR -DPROD (PR,UAEF(1 >)*PXR 
;>YR  =  OPROD(PR*UAEF(2)  >+PYR 
PXI=OPROO(PI*  UAEF<1  )  >*PXI 
PY  I  =DPROO  (PI*UAEF(2))>PYI 
500  CONTINUE 
525  CONTINUE 

IF  (YC.ECUO  .0),  GO  TC  200 
IF  (YC.GT.C.O)  THEN 
YC=-YC 
GO  TO  TOO 
ENOIF 
GO  TO  200 

550  IF  (XC.GT.C.O)  GO  TO  150 
600  CONTINUE 

C  FAR  FIELO  ELECTRIC  FIELD  (SQUARED)  CALCULATED 

EVE.RT2=(PYR**2*PYI**2)*CTHETA**2 
EHOR2=(PXR**2+PXl**2)*CTHETA**2 
C  PATTERN  CALCULATED  -  NORMALIZED  TO  VALUE  AT  THETA=?HI=0 

PVERT=0B(EVERT2/PMAX) 

?HOP=OB(EHOR2/PMAX> 

C  VARIABLES  OF  SUMMATION  OUTPUT  FOR  USE  IN  PROGRAM  MAIN2 

WRITE (12*4800  )THET AD» PXR »PX I  * PYR » PY I 
♦800  FORMAT  (F7 . 3* / » ♦( E ♦ 0 .33 »/ ) ) 
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non 


C  OUTPUT  TO  TERMINAL  -  MA v  HE  DELETEED  IF  NOT  OESIRED 

WRITE  < *»5OO0 )TH£TAO»PVERT»PHOR  »DSHIFT ,PXR  »PXI »PY  R  »PY  I 
C  ANTENNA  PATTERN  OUTPUT  TO  FILE  PATTERN 

WR I TE ( 10  »  5C  0 0 )  THETAD»PVERT  »PH0R 
5000  FORMAT ( 7< F7 . 3 *2X ) ♦ / • 4E1 0 . 3 ) 

650  CONTINUE 
703  CONTINUE 
STOP 
END 

THIS  FUNCTION  ALLOWS  ANY  PARABOL ICALLY  TAPERED  APERTURE  FIELD 
TO  3E  MOOELEO  WITH  PARAMETERS  IN  ANTENNA  FILE.  MAY  BE  REPLACED 
WITH  ANY  OESIRED  ILLUMINATION 
FUNCTION  EF  IELD (R *R  A  *C  » N) 

IF  (N.EQ.0.0)  THEN 
EFI ELD=1. 0 
GO  TO  10 
ENOIF 

£FIELD=C* (1-C)*(1-<R/RA)**2)**N 
10  RETURN 

END 

FUNCTION  DB ( X ) 

Y=ARS(XI 

IF  ( V.GT.1.0E-O8)  GC  TO  10 
0B=-80 .0 
GO  TO  20 

10  DB=10.O*ALOGlOf Y) 

20  RETURN 
END 


P  R  0  5R  A  M  MAIN2 

DIMENSION  U!7>  tOX!7)tDY!  7)  *U  AEF  <2  )  ♦  X  D  <  10  )  ♦  TD  ( 1 3  >  »  3P  ( 1 0  > 
♦♦0SHlFT!4)f0ISCR!4) 

REAL  MN!4> 

33U3LE  PRECISION  PXR f PXl tPYR fPY I 
OPEN! 7, FILE=* ANTENNA* ) 

3PEN!8fFILE=*  ANGLES  * ) 

3PEN!9,FILE=*  DISCS* > 

3PE *H  1 3  »FILE= ‘PATTERN •  ) 

3PEN(11#FILE=*PA* ) 

DATA  DX(1)*DY(1)*DY(2)  t  DY! 3 ) /4*  0 .0/ 

REMIND! 7) 

REW I  NO (8 ) 

REW  I  NO  (  9) 

REMINO(ll) 

READC7,10  00)SIDE»R A DIUS » FREQ »F » BLOCK »C»N»PMAX,  1 01  SC 
10  00  FORMAT! /» S( /, El 4. 7« / ) « /I2 1/ . /♦ E 14 . 7 » / ,/ .  13 > 

READ! 8t200C  >PHI t  NANG 
2000  FOR  MAT !/f/*F10»7»/# /»I3»/) 

30  50  IR= 1» IDISC 

READ!9f2600)XD!IR)tYD!IR) tDISCR !IR) t DSHIFT ! I R ) 

2600  FORMAT! 3!/,E14.7t /) »/fE14.7) 

DEL  T  A=0SHIFT !IR >*3. 0E+O8/FREQ 
RP=SQRT!XO!IR>**2*YD!IR>**2> 

MN! IR)--RP/!2.0*F) 

3P!IR)=!RP**2*4.0*F**2I/!2.0*RP>- 
♦DELTA*t2.0*F-MN!IR>*RPJ/!RP*SQRT!  1 . 0*>MN!  IR  > **2)  ) 

50  CONTINUE 

PIE=3. 14159 26535898 

9ET  A=2.G*PIE*FREQ/3.0E+G8 

PHI  =PHI*P  IE/180  .0 

M!l)  =  258. 0/1003  .0 

W!2)=125. 0/1008.0 

U! 3)=W!2) 

U(4)=W(2) 

W!5)=M!2) 

M ! 6 ) =W ! 2 ) 

W(7)=U(2) 

OX! 2>=SIDE*SQRT!14.0>/5.0 
OX!3)  =  -OX!2) 

DX! 3>=SIDE*SQRT! 14. 31/10.0 
OX!7>=0*! 5) 

OX ! 4)=-0X!5> 

DX!6)=0X!4) 

3Y!  4 )  =  SIOE*SQRT! 42. 01/10.0 
DY ! 5)=0Y! 4) 

DY  !  ->)s-DY  !4  ) 

0Y!7)=0Y!S) 

CPHI=COS!PWI> 

SPHIsSIN!PHI) 


yincr=3.o*s  ide 
<  INCR=SQRT(  3.  G)  *SIDE 
N1=INT(RA0IUS/XINCR > 

DO  700  L=1  ♦ NANG 

RE AD<  11  *300 G)THET AD  tPXR  »PXI »PYR  *PYI 
FORMAT(F7.3*/ »3(E 40,3 3,/) »E4 0.331 
THETA=THETA0*PIE/130.Q 
CTHETA=COS(THETA) 

STHET  A=SI N  <  THET  A) 

XC=XINCR/2.0 
00  600  IX=1 »2*N1*1 
XC  =  -XC*XINCR/2.0 
I TEST=MOD (I  X*  2) 

IF  (ITEST.GT.O)  THEN 

rc=o.o 
rcc=o .o 

SO  TO  300 
ELSE 

YC=YINCR/2.0 
YCC=YC 
GO  TO  300 
END  IF 

YC  =  YC*YINCR 
GO  TO  300 
XC=-XC 

rc-rcc 

GO  TO  300 
YC=-YOYINCR 
CONTINUE 
I DT  EST-0 

DO  350  ID=1 tIDISC 

DIST  =  SQRT  (  (XC-XDCIO) )**2+< YC-YD( ID) )**2 

IFIDIST.LT. DISCRilD))  GO  TO  375 

CONTINUE 

IDT£ST=1 

CONTINUE 

RPRIME=SQRT(XC*XC*YC*YC) 

IF  (RPRIME.LT. BLOCK)  GO  TO  100 
IF  (RPRIME.GT .RADIUS)  GO  TO  550 
IF  <IDTEST.C7T.O)GO  TO  525 
DO  500  1=1*7 
X  =  XC+DX  <1 ) 

Y=YC+DY(I ) 

RPRIME=SQRT  <X*X>Y*Y ) 

PHI P= A  T  AN2 ( Y  »  X ) 

THETAP=2.0*ATAN(RPR IME/(2.0*F) ) 
CPP=COS(PHIP) 

SPP=SIN(PHIP) 

CTP=COS(THETAP) 

STP=SIN(TH£TAP) 


0=S1RT(1.0-SPP*SPP*STP*STP) 

UAEF(1 )=(3PP*CPP*<1 .0-CTP>/D) 

UAEF (2  )  =  (SPP*SPP*CTP*CPP*CPP) /D 
C2=EFI  ELD  (R PR IM E# RADIUS «C»N)*U<  I ) 

ARG=SET  A* RPRI ME*S THET  A*COS< PHI-PHI? ) 

XP=(RPRIME-fiP(ID>  >*MN(ID> 

3  =  SQRT(XP**2«-RPRIME**2 >-*XP 
SHIFT=(2.0*F-S)*8ETA 
ARGO=ARS- SHIFT 
PR=C2*C0S (ARGDI 
3I=C2*SIN(ARGD> 

PR0=C2*C0S( ARG) 

31 0=C  2*SI N( AR  G> 

3XR=DPR0DCPR»UAEF(1  )  )+PXR-DPROD (PRO »U AEF ( 1 ) > 

PYR=DPROO  (PR»UAEF<  2 ) > ♦PYR-DPROD (PRO  »UAEFC2> ) 

3  XI  =  OPR 00 (PIfUAEF(l))+PXl-DPROO(PIOtUAEF(l)) 

3YI=DPR 00(PIt  UAEF(2  > >+PYI-OPROD(PIO  *UAEF(2) ) 

500  CONTINUE 
525  CONTINUE 

IF  (YC.EQ.O.O)  GO  TO  200 

IF  (YC.GT.0.01  THEN 

YC  =  -YC 

GO  TO  300 

ENOIF 

GO  TO  200 

550  IF  (XC.GT.O.OI  GO  TO  150 
600  CONTINUE 

EVERT2=<PYR**2»PYI**2)*CTHETA**2 
EH0R2=(PXR**2+PYI**2)*CTHETA**2 
PVERT=0B(EVERT2/PMAXI 
PH0R=DB(EH0R2/PMAX  > 

OUTPUT  TO  TERMINAL  -  MAY  BE  DELETED  IF  DESIRED 
URITE(*»500ClTHETAD*PVERT»PHOR»DSHIFT,PXR,PXI*PYRfPri 

ANTENNA  PATTERN  OUTPUT  TO  FILE  PATTERN 
UR  I  TEC  10 *53 001  THET AD»PVERT ,PHOR 
5000  FORMAT! 7(F7.3»2X) ,/ |4E10. 3) 

650  CONTINUE 
700  CONTINUE 
STOP 
ENO 

FUNCTION  EFIELD(RtRAtCiN) 

IF  (N.EQ.O.CJ  THEN 
EFIELO=1.0 
GO  TO  10 
ENOIF 

EFIELD=C* (1-C)*(1-CR/RA)**2)**N 
10  RETURN 
END 

FUNCTION  D8(X ) 

T=ABS(X) 


FILE^ANTENNA 

SIDE  (HEXAGON  RADIUS) 

.01 

RADIUS  (RADIUS  OF  THE  APERTURE  PLANE) 

0.6096 

FREQ  (OPERATING  FREQUENCY) 

2. 8E*09 

F  (FOCAL  LENGTH  OF  DISH) 

•  6858 

BLOCK  (RADIUS  OF  AREA  BLOCKED  BY  FEED) 

.0762 

C  (FEED  PARAMETER-ILLUMINATION  AT  EDGE  OF  DISH) 

0.1 

N  (FEED  PARAMETER-POWER  OF  PARABOLIC  TAPER) 

1 

PMAX  (NORMALIZATION  VALUE  FOR  MAINOS) 

5. 736891  IE +  06 

IDISC  (NUMBER  OF  OISCS  ON  DISH) 


XI 

-0.2413 

ri 

o.o 

DISC  RADIUS  (DISC*) 

0.08S9 

DFUNC1 

0.2625 

X2 

0.2413 

Y  2 

0.0 

DISC  RADIUS  (DISCR) 
0.0889 
DFUNC2 
0.260 
X3 

-0.1143 

Y3 

0.0 

DISC  RADIUS  (DISCR) 

0.03175 

DFUNC3 

0.306 

X  4 

0.1143 

Y  4 

0.0 

DISC  RADIUS  (DISCR) 
0.03175 
DFUNC4 
0.201 


FIL£= ANGLES 

PHI  (ANGLE  THE  PATTERN  IE  BEING  TAKEN  ACROSS) 

00.0 

NANG  (NUMBER  OF  ANGLES  THAT  POLL OW-< THET A )  ) 

4 

THETA  (ANGLE  THAT  THE  PATTERN  IS  BEING  COMPUTED  FOR  IN  DEGREES) 
10.0 
11.0 
12.0 
13.0 


V 
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EXAMPLE  OF  FILE  PA  AS  GENERATED  BY  MAIN1T.  AND  USED  BY  MAINT2 

0.0  0 

•3331394 33 2 3 391494 79 59 47 38 018 94 170E-14 

. 2 39518 0 78 471 87817798 7388243795982E+0 4 

• 

1.00 

.  30 3 43 44 0 32 32 0 6 79 69 C 57 7 0207621 6 272E-1 4 
. 1319356097955743833732682221292 49 E- 14 
.2309 1067709729478390 21 5621 4516955E *04 
. 24  72  055 59 10 29824 09781 9 99 9558 00 99 IE- 12 
2.00 

•2262548694308258 48 35114493 44 3010 4E -14 
.23710637427334 42635793 392 17 79371 7E- 14 
.20 643 1893 051 3714 662686 90 54 7520 6 09E+0 4 
.56672 01 3588103 790868851 4329420 500E-12 
3.00 

.10 938 9632 563 05 4368 53 38 24272258 145E-1 4 
.295311012970 18 8786 5468 83 388580 703E- 14 
. 1698 3 1 78 6582 48 56 74 670 7700 7CO 38 8 04E+0 4 
.6393110997214251 9850 6494369047627E-12 
4.00 

1526  0  9  307  2352721976640  618320  455  21E-15 
297665388428 1692782618 92702558514E-14 
12b4 95331 863592154 7468 4460 9862255E+0 4 
61 1266 318130 57532991 18 5853C351 342E-12 
5.0  0 

1233 76 329 65914 9394 4569 839521 42920E -14 
24 R279 5704 7449 6413 1149 5879 7189 792E- 14 
82 366831 08 12 485882 62 35 35 3205 054 15E+0 3 
27 51 0 5 292529 91 48 359964 66 28 320453 8E -12 
6.0  0 

194441 439 3 7694 7767667 48 00 70281 599E- 14 
1625  0  3988  0  5629  52  413228  726  677838 40E-1 4 
42825 7650 328 35 39 2918 6447 8 143 53 51 IE ♦03 
13332 7738 5R347 67 75350 1746 594 127 05E- 12 
7.0  0 

21 366 0 377255893 12770 7 3 46 31 25 18 29 8E- 14 
6233 1 1621 12040 32830 02431 5675633 48E- 15 
1176834013 98 79 7257294430 7821550 40E* 03 
60 60 79 0 06860 94 20 25 7477 9 748 12 21 498E-1 2 
3.00 

l',79189359  75  773719771S70339355  368E-14 
29293011922909C187282609261928311E-15 
8912713399  06  01263596  40  710  93571 50 2E* 0  2 
10 979 12895 7781 5938 1931 45524549 3 19E- 11 
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